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Transport  Properties  of  Dense  Laser  Plasmas 
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2 Fachbereidi  Pliysik,  Universitat  Rostock,  Universitatsplatz  3,  18051  Rostock,  Germany 

Generalized  kinetic  equations  are  used  to  investigate  transport  properties  of  dense  plasmas 
in  strong  laser  fields.  Quantum  statistical  expressions  are  derived  for,  e.g.,  the  collision 
frequency  to  account  for  nonlinear  field  phenomena  as  well  as  strong  correlation  effects.  The 
influence  of  these  effects  on  collisional  absorption  is  discussed. 


1.  Introduction 

The  impressive  progress  in  laser  technology  makes 
femtosecond  laser  pulses  of  very  high  intensity  avail- 
able in  laboratory  experiments.  If  solid  targets  are 
irradiated  by  such  laser  pulses,  dense  plasmas  can  be 
created  relevant  for  astrophysics  and  inertial  confine- 
ment fusion.  To  understand  transport  properties  of 
intense  laser-plasma  interaction,  a theoretical  mod- 
elling of  dense  nonideal  plasmas  in  strong  laser  fields 
is  needed. 

One  of  the  important  mechanisms  of  energy  deposi- 
tion is  inverse  bremsstrahlung,  i.e.,  laser  light  absorp- 
tion via  collisional  processes  between  the  plasma  par- 
ticles. In  strongly  ionized  plasmas,  this  absorption 
process  is  essentially  governed  by  the  electron-ion  in- 
teraction usually  described  in  terms  of  the  electron- 
ion  collision  frequency. 

A lot  of  work  has  been  done  up  to  now  to  consider  the 
electron-ion  collision  frequeny  and  the  dynamic  con- 
ductivity, respectively,  for  laser  plasmas  under  dif- 
ferent conditions.  The  well-known  classical  theories 
were  developed  in,  e.g.,  [1,  2,  3].  Basic  equations  of  a 
rigorous  quantum  kinetic  approach  to  laser  plasmas 
were  derived  in  Ref.  [4]. 

We  consider  a fully  ionized  plasma  under  the 
influence  of  a spatially  homogeneous  electric  field. 
Starting  point  is  a generalized  kinetic  equation  for 
the  gauge  invariant  Wigner  distribution  function 
which  follows  from  the  time  diagonal  Kadanoff- 
Baym  equation.  The  central  quantity  in  the 
collision  term  is  the  self  energy  function.  Powerful 
schemes  are  available  to  determine  appropriate 
approximations  for  the  self-energy  function  taking 
into  account  nonlinear  field  dependence  as  well  as 
many-body  and  quantum  effects  relevant  for  high 
density  plasmas. 

2.  Collisional  absorption 

As  we  are  interested  here  in  the  collisional  absorption 
by  the  plasma,  it  is  obvious  to  start  from  the  balance 
equation  for  the  energy  and  for  the  electrical  current 
resulting  from  the  kinetic  equation.  As  the  collison 
integral  is  a non-Markovian  one,  the  energy  balance 
reads 

dWk'n  dWpat  . _ 

-dT  + ~ir  ~J'E’  (1) 

i.e.,  the  change  of  the  total  energy  of  the  plasma 
particles  is  equal  to  j • E that  is  in  turn  the  energy 


loss  of  the  electromagnetic  field  due  to  Poynting’s 
theorem. 

For  the  calculation  of  the  collisional  absorption,  we 
start  from  the  general  balance  equation  for  the  cur- 
rent density  of  species  a.  It  follows 


dt' 
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where  Vab  ( q ) is  the  Fourier  transform  of  the  Coulomb 
potential.  The  collision  term  is  now  expressed  by  the 
correlation  function  of  the  density  fluctuations  given 
by  ihL<b(t,  t')  = (5pb{t')5pa(t)).  For  this  function  an 
appropriate  approximation  has  to  be  found.  In  Born 
approximation  the  current  can  be  expressed  by  the 
Lindhard  dielectric  function.  For  a harmonic  field 
E = Eocoswt,  the  field  dependence  is  described  by 
a Fourier  series  in  terms  of  Bessel  functions.  This 
scheme  allows  to  investigate  the  current  density  and 
the  energy  dissipation  in  dense  weakly  nonideal  laser 
plasmas  including  dynamic  screening,  quantum  ef- 
fects as  well  as  nonlinear  field  phenomena  such  as 
higher  hamonics  and  multiphoton  processes.  The  de- 
tails of  the  investigation  are  presented  in  Ref.  [6]. 

To  include  strong  correlations  one  has  to  go  beyond 
the  Born  approximation.  An  approximation  which 
accounts  for  strong  correlations  in  the  ion  and  elec- 
tron subsystems,  we  derived  in  Ref.  [7].  Now,  the 
resulting  expression  for  the  current  is  given  in  terms 
of  dynamical  structure  factors  and  the  exact  den- 
sity response  functions.  Again,  the  field  dependence 
is  described  by  a Bessel  function  expansion  which 
causes  the  nonlinear  field  effects  mentioned  above. 
An  important  quantity  is  the  cycle  averaged  dissipa- 
tion of  energy  (j  • E).  In  the  case  of  high-frequency 
fields,  we  get 
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with  ee1{q,;to)  = 1 + (h2  e2e)  / (eQq2)  £fe(q;w)  being 
the  dielectric  function  with  the  density  reponse  func- 
tion of  the  electron  gas  £ee(q).  As  high  frequency 
fields  are  considered  the  static  ion-ion  structure  fac- 
tor Si,(q)  appears  in  (3)  instead  of  the  dynamical 
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one.  Furthermore,  Jn  is  the  Bessel  function  of  nth 
order  and  vq  = eEo/meco  is  the  so-called  quiver  ve- 
locity. 

Often  the  electron-ion  collision  frequency  vei  is  dis- 
cussed which  is  defined  for  the  high-frequency  case 
by  (wp  - plasma  frequency) 
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Using  the  Lindhard  dielectric  function  and  ne- 
glecting ion-ion  correlations  (Sa=l),  the  results 
follow  we  derived  in  Ref.  [6].  Those  results  have 
a similar  form  as  that  of  the  nonlinear  Dawson- 
Oberman  model  [1].  However,  in  our  formula,  the 
dielectric  function  is  given  by  the  quantum  Lind- 
hard form,  whereas  the  dielectric  theory  of  Decker 
et  al.  leads  to  the  classical  Vlasov  dielectric  function. 


3.  Numerical  results 

We  compare  the  results  obtained  in  the  frame  work 
of  the  quantum  kinetic  approach  with  the  results 
of  the  classical  theories.  In  Fig.  1 the  collision 
frequency  is  shown  as  a function  of  the  coupling  pa- 
rameter T = (e2 / 47re0) / dkDT  with  d = (47rn/3)-1/3 
for  a strong  field  ( vo/vtn  > 1)-  Results  using  the 
Lindhard  dielectric  function  and  Su( q)  = 1 are 
given  by  the  upper  curve.  In  the  case  of  static 
screening  the  collision  frequency  is  reduced. 


of  the  coupling  parameter  F for  a hydrogen  plasma  in  a 
strong  laser  field.  Comparison  is  given  with  the  theory  of 
Decker  et  al.  and  with  the  asymptotic  formula  of  Silin. 

For  small  values  of  T (weakly  nonideal  plasma),  the 
classical  dielectric  theory  of  Decker  et  al.  [1]  and  the 
classical  asymptotic  formula  by  Silin  [2j  give  almost 
the  same  results.  However,  with  increasing  T,  the 
classical  curves  show  a sharp  drop  down.  This  be- 
haviour results  from  a cutoff  procedure  at  large  wave 
numbers  used  in  the  classical  theories.  Such  a cut- 
off is  avoided  in  our  approach  automatically,  and  the 
range  of  applicability  is  extended  to  higher  values  of 
the  coupling  parameter. 

As  discussed  above,  the  expression  (3)  generalizes 
the  results  to  dense  quantum  plasmas  including 
strong  correlation  effects. 


of  the  coupling  parameter  T.  Ion  structure  factors  in 
different  approximations:  HNC  and  Debye. 

Let  us  consider  the  influence  of  the  ion-ion  corre- 
lations. In  Fig.  2 the  collision  frequency  is  shown 
using  different  approximations  for  the  static  ion-ion 
struct, nr  factor  Su  (in  the  dielectric  function,  we  in- 
cluded local  field  corrections  [7]).  Here  Vo/vth  = 0.2, 
i.e.  a smaller  field  strength  is  considered.  Inclusion 
of  the  structure  factor  in  hypernetted  chain  (HNC) 
approximation  decreases  the  collision  frequency  for 
small  and  moderate  coupling  up  to  T « 5 whereas 
an  increasing  behaviour  follows  for  high  values  of  the 
coupling  parameter.  As  expected,  the  Debye  approx- 
imation can  be  applied  only  for  weak  coupling.  Good 
agreement  with  data  obtained  by  simulations  can  be 
found  in  the  range  T < 1 [6). 
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